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We consider a class of modified gravity models where the terms added to the standard Einstein- 
Hilbert Lagrangian are just a function of the metric only. For linearized perturbations around an 
isotropic space-time, this class of models is entirely specified by a rank 4 tensor that encodes possibly 
time-dependent masses for the gravitons. This tensor has the same symmetries as an elasticity 
tensor, suggesting an interpretation of massive gravity as an effective rigidity of space-time. If we 
choose a form for this tensor which is compatible with the symmetries of FRW and enforce full 
reparameterization invariance, then the only theory possible is a cosmological constant. However, in 
the case where the theory is only time translation invariant, the ghost-free massive gravity theory is 
equivalent to the elastic dark energy scenario with the extra Lorentz violating vector giving rise to 
2 transverse and 1 longitudinal degrees of freedom, whereas when one demands spatial translation 
invariance one is left with scalar field theory with a non-standard kinetic term. 



I. INTRODUCTION 



The realization that the Universe appears to be accelerating has fuelled the search to alternative theories of gravity 
PQ as a possible explanation for what has become called the dark sector. In this paper we will focus on what is the 
simplest subset of such theories, in which the dark sector Lagrangian is only a function of the metric (and no extra 
derivatives thereof) following the approach discussed in [2J. The action for this type of theory is given by 



S 



R + WnGC m -2£ i {g llv ) 



(1.1) 



If T M „ is the energy-momentum tensor for the matter sector Lagrangian £ m and is that associated with the 
dark sector Lagrangian £ d , then the Einstein equation is — SttGT^ + U^. Typically, we will be interested 
in spacetimes which are isotropic where — pu^u v + Pj^ can specified in terms of a density, p, and pressure, 



P = wp. There arc two classes of theories which can be described by ( 1.1 1: clastic dark energy and theories of massive 
gravitons. Typically, theories of massive gravitons have been considered as a fundamental theory around Minkowski 
space-time, but it is also possible to think about masses for the gravitons as being induced by a some unknown 
effective physics encoded by C d . 

The study of massive gravity theories (see, for example, [3111]) has a long history. This started with the linearized 
theories of Pauli and Feirz [5J, progressing to the studies of Boulware and Deser [6J. It has received a new lease of 
life in recent times with the proposal of non-linear dRGT massive gravity theory [7lll0| and its connections to the 
Vainshtein screening mechanism |11H14| . Massive gravity theories are built upon the pretext that the resulting theory 
should be "ghost-free" IT5rf2"l"] , and have begun to be studied in cosmological backgrounds [2"2"rf3"4"] . Such theories 
are usually presumed to be Lorentz invariant which leads to the Pauli-Fierz tuning. Giving up Lorentz-invariance is 
another way to remove ghosts from the theory, as pointed out by (35] [36] and further studied in (3[ I37H40] . 

Elastic dark energy (EDE) is an idea which has been developed from relativistic elasticity theory [4"TH4"T] . The 
basic concept is that the stress-energy component responsible for the dark energy has rigidity which stabililizes 
perturbations that would, if modelled as those of a perfect fluid, give rise to exponential growth in the density 
contrast. The framework was adapted for cosmological purposes in [48, 49 in order to provide a phenomenological 
model for domain wall as an explanation for accelerated expansion which have P/p — w dv = —2/3. However, in 
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principle the equation of state parameter, w, is allowed to take "any" value so long as the rigidity modulus, fj,, 
is sufficiently large. Indeed, the theory is well-defined in the limit w — > — 1 where the elastic medium becomes a 
"cosmological constant" and w — > and fi — > which corresponds to cold dark matter. The standard assumption, 
which we will use in this paper, is that the elasticity tensor is isotropic, but one can also construct anisotropic models 

EDH5U. 

The aim of this paper is to point out the connections between linearized massive gravity theories, EDE and the 
framework for linearized perturbations in the generalized models for the dark sector discussed in [2, 52]. The reason 
for this connection is that, at linearized order, one can represent all possible Lagrangians by a generalized function 
which is quadratic in the fields. In the specific case we are concerned with here, this is just a quadratic function of the 
metric which is parameterized by a rank-4 tensor. This has the same symmetries as an elasticity tensor, suggesting 
an interpretation of massive gravitons as creating an effective rigidity of spacetime. This tensor can split in a way 
which is compatible with the symmetries of the FRW spacetime, that is more general than the usual Pauli-Fierz case. 
We will make a survey of the possible mechanisms by which ghost modes can be removed, both Lorentz invariant 
and violating. The ghosts are associated with a breaking of reparameterization invariance and we show that its re- 
imposition leads to 3 interesting sub-cases, one of which is compatible with Lorentz invariance which is a cosmological 
constant and the other two which violate either time or spatial translation invariance. The one which violates spatial 
translation invariance happens to be the EDE model which we see is a Lorentz-violating ghost-free massive gravity 
theory. 



II. THE GENERAL "METRIC ONLY" THEORY 



The action which will give linearized field equations for the perturbed field variables is given by 



S {2 } 



j4 / ' 

d x^J-g 



2 i?+ 167rG0 2 A„ - 2£ {2 



(2.1) 



We use 2 Q to denote the second measure-weighted variation of the quantity Q, defined as 2 Q = v /=^ 2 (v / ~ff Q)- 
C {2 ) is the Lagrangian for dark sector perturbations, given by 



A 2 > = eW^Sg^Sgap, 



(2.2) 



where Sg^ lL , is the metric fluctuation. This clearly looks like a mass term for the metric perturbations. The tensor 



(2.3) 



is the mass matrix determining how the components of Sg^ mix to provide the mass; all linearized massive gravities 
are encoded by choices of W. Hence, the complete linearized theory we study is 



S m -J 



j4 / ' 

d xy/-g 



2 R + 16vrG0 2 /: m - \W^8 gtlu 5g a0 



(2.4) 



This theory contains metric fluctuations which have a kinetic term, and a mass term. We will see that the spatial 
components of W can be interpreted as an elasticity tensor. 

To isolate the degrees of freedom in the theory, Sg^ can be decomposed as 



(2.5) 



In the parlance of [42, 47, 53, 54 h^ v is the Eulcrian metric perturbation and is a vector field representing possible 
coordinate transformations. In standard General Relativity, the action is independent of but in more general 
theories this can become a physical field. This formulation is equivalent to what is sometimes called the Stuckelberg 
trick 4, 18, 55]. Inserting (2.5) into the action (2.4) and integrating by parts reveals that 



S 



(2) 



^4 

CI X\ 



(S E G^ - 8ttGS e T^ - S E U^)h^ + 2£ (A1 <S E (V„)[/n 



(2.6) 



where the variational operator "5 E " denotes that the quantity is evaluated with the metric perturbation variable 
(rather than Sg^), and where we have defined the perturbed dark energy-momentum tensor 



5 E U^ 



(2.7) 
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It is now a simple matter to obtain the functional derivatives of the action with respect to the perturbed metric 
and ^-fields, 

S {2} = S E G^ - 8ttGS e T^ - 5 E U^ = 0, (2.8a) 

J^S {2} = 8 E {V V U^) = 0. (2.8b) 
The variational principle was used to demand that these expressions vanish, yielding the perturbed gravitational field 



equations and perturbed conservation equation respectively. Using (2.7) to evaluate (2.8b I yields 

L^V^^p + (V M W^)V Q ^ + Flaunt? = -|((V M W^)^ + P^V^Kp), (2.9) 
where we defined the effective metric L^ ua @ and derivative-coupling pv va P terms, 

jjiua.fi — _|_ jjpv gUp _ 2JJot{^g v )fi pfJ.va/3 = yyfiuafi _|_ jjafi gfiu _ 2g v ^\J a ^ (2.10) 

We impose spatial isotropy upon the background with the (3+1) decomposition and in doing so we will obtain the 
most general linearized massive gravity Lagrangian compatible with spatial isotropy of the background. We foliate 
the 4D spacetime by 3D sheets with a time-like unit vector, u^, being everywhere orthogonal to the sheets. The 
4D spacetime has metric <? M „, and the 3D sheets have metric 7^. The (3+1) decomposition of the 4D metric is 
9fii/ = 7fiv — u^u v , where u^u^ — —1, 7 A " y u A1 = 0. This structure provides an extrinsic curvature tensor = 
on the 3D sheets, given by K^ v = V M u„ and satisfying u^K^ v = (the extrinsic curvature tensor is given by 
K^u — \K"tnv)- We define "time" and "space" differentiation as the derivative operator projected along the time and 
space directions, 

= u"V^, = 7%V„V- (2-11) 

Using this technology we decompose the gradient of a scalar into two orthogonal terms, 

V M ^ = -^u M + V M V- (2-12) 
This enables us to find the values of two useful "kinetic scalars" , 

Wi/jV^ = -ip 2 + V^V^V, (2.13a) 

Dtp = V^V^V = -$ + V^V^V- (2.13b) 

The last term of each expression simply selects out the spatial derivatives of the scalar field. Another useful application 
of the (3+1) decomposition is to find all the freedom in a tensor which is compatible with spatial isotropy of the 
background spacetime. 

We use the (3+1) decomposition to isolate the components of the perturbed metrics by writing 

Sg^ = 2<f> Utl u u + 2N iM u u) + H af}1 a (2.14a) 

hpv = tyu^ + 2n^u v ) + h afj j a ^ v , (2.14b) 

and we isolate the time-like and space-like components of the vector field via 

£ M = -Xu^ + u}^, (2.14c) 

where N^u^ = n^u^ — 0, u^H^ — uNi^ — and u^uj^ = 0. 

In [2 j we showed that the general decomposition of the mass-matrix W compatible with spatial isotropy is 

yynvafi _ AyvU^U^u"^ + B W (u^u"^ + U a U P 'f U ) + 2C W ('f^U^vP + "f'^U^u") 

+D w Y v l afi + 2E w Y (a l P) \ (2-15) 
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where there are only 5 free functions which only depend on time. Using the mass matrix (2.15) and (2.14a) in the 
Lagrangian (2.2) yields 



8£ {2} = 4A w <S> 2 +4B w <$>H + 2C w N a N a +D w H 2 + 2E w H al3 H al3 . 
This can be written in terms of "graviton masses" (see e.g., [3]) where the free functions {^4yv>! 

2£ {2} = Tnl(Sgoa) 2 + 2m\(5 gm) 2 - m%(8g i j) 2 + m 2 (5g ti ) 2 - 2m\5g m 5g il 



(2.16) 

. , Eyv} are given by 
(2.17) 



with A 



w 



'()• 



D 



w 



-2ml O 



vv 



2 D w — 4m 2 and Ew = —2m\ 
, Ey\>} we are actually talking about the graviton masses m 2 , albeit ones that 



2m\. One should keep in mind, therefore, 



that when we talk about the {Aw 
depend on time. The values of these masses for a "Goldstone" theory are given in [56] and those which are induced 
by perturbations in scalar fields in [2J. 
Using (2.14) to evaluate (2.5) yields 



•X) 



(2.18) 



where w a = (&> a — ^Kuj a ). Substituting (2.14a) into the action (2.4) one finds the absence of 4> 2 and h 2 a terms in 
the kinetic part of the Einstcin-Hilbcrt Lagrangian; this can also be seen in results given by [31 [SHI EH] and in the 
ADM formulation 57J. <j) and n a are now Lagrange multipliers whose equations of motion are constraint equations, 
allowing them to be eliminated. Using re-definitions of the coefficients, we can effectively set (f> — and n a — which 
is equivalent to choice of the synchronous gauge. We will make this choice in what follows. 

The two independent components of the equation of motion (2.9), after inserting the (3+l)-decomposition, are 
given by 

[A w + p] x + [Aw + n(4A w +p- 3P)] ± + [P + Cw] V 2 X 

- [U{3P + 2p-A w + 3B W ) - (2A W + 5p + 3P- 3B W - 9D W - 6E W )H 2 

+ (2p + 3B w -A w )^+p] x 

- [B w + H(3B W + 3D W + 2E W - 2P)} d l uj i + [B w + C w ] diLu 1 



= ^[Bw 



H(3B W + 3D W + 2E W - 2P)] h+- [B w - P] h, 



(2.19a) 



[P- 



C w p l - [Cw+U^Cw- p + 3P)]Co l - [E W -P]V 2 J- [E w + D w ]d l d k uj k 
+ [C w - H(3D W -B w + 2E W - 4C W - 2P)} d\ + [B w + C w ] d l X 

= -[P- Ew] djtii + \[P + Dw] d l h. 



(2.19b) 



p and P are the density and pressure coming from the dark fluid (i.e. the components of the background dark energy 
momentum tensor U^). The benefit of using the (3+1) decomposition has become apparent: we are able to identify 
the degrees of freedom. There are the tensor degrees of freedom h+ and h x which are present in standard General 
Relativity, a vector degree of freedom u/, that can be split into a longitudinal (scalar) and two transverse (vector) 
degrees of freedom, and a scalar degree of freedom x- Therefore, prima facie there are 6 extra degrees of freedom. As 
we will discuss below either x or w * can be a ghost and therefore the coefficients must be chosen to suppress one or 
both of them. In the case where x is the ghost then there are 5 degrees of freedom with those in oj 1 being split into a 
longitudinal (scalar) and two transverse (vector) degrees of freedom. If oo l is the ghost then there are only 3 degrees 
of freedom. 



III. MECHANISMS FOR THE ELIMINATION OF GHOSTS 



From (2.6) we see that the kinetic terms of the w M and x fields enter the theory via 
1 



(A w + p)x! + (Cw + P)V AlX V^x 



S {2} D / d xy/^g 1 



+(C W - p)^^ + (D w + P)(V^) 2 + 2(E W - P)V^VV 



(3.1) 
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Let us now focus on the standard scenario of perturbations around Minkowski space-time when both p = P = 0. If 
Ay\> > 0, then Cyy < is required for \ to have a kinetic term with the "proper sign". But if this is the case then w M 
has a kinetic term with the "wrong sign" (the same is true if Ayj < 0). Hence, one of %, w M must be a ghost. There 
are a few ways to get out of this. 

First, one can make the coefficient of x 2 vanish by setting Ayy = 0, which removes \ as a propagating mode and 
the equation of motion is a constraint that can be enforced in a Lorentz invariant theory. Secondly, one could set 
Cyy = 0, since that would remove as a propagating mode. Finally, one could set x = directly which requires a 
breaking of Lorentz invariance (since we will be manually forcing one of the four components of a 4- vector to zero). 

When Ayj = 0, there is no \ 2 term in the Lagrangian and the equation of motion simply becomes a constraint 
equation specifying the value of x from the other field variables. This can be back-substituted into the action so 
that the theory explicitly does not contain the %-field. From our presentation it is clear that in this case the ghost 
can be identified with the time-like degree of freedom x- We have been able to deduce this since we used a (3+1) 
decomposition. Performing a trans verse- longitudinal decomposition does not aid the identification of the ghost. 

If we choose the 5 parameters in (2.15 1 to be given by Aw = X + Y, £?yy = —X, Cw = — \Y , Dyj = X, Eyj = \Y , 
then the theory is Lorentz invariant and the mass-matrix is be given by 



W' 1 



xgnr p 



Yg^ a g 



(3.2) 



where X, Y are two parameters that are dependent on background field variables only. The standard route for isolating 
the ghost in the Lorentz invariant theory [55] decomposes £ M into its transverse and longitudinal modes as 



C> = G< + v^k, 



(3.3) 



where VmC^ 1 = and n is a scalar field. For Minkowski background spacetime, inserting (2.5), (3.3) and the mass- 



matrix (3.2) into the Lagrangian (2.2), whilst assuming that X,Y are constants, yields 



C {2} = l -{Xh 2 + Yh^h^) + l -Yh^d { ^ v) + l -{XhU K + YhTd^duK) + \Yd»Cd {p Q v) 



+Yd^C ) d fl d vK +\{X + Y){Unf 



(3.4) 



This expression has made the ghost problem manifest in the Lorentz invariant language. The existence of the last 
term, (X +Y)(Dk) 2 , means that ghosts are inevitable (see e.g. 15, 58, 59]). The cure is to set X = —Y, which removes 
the problematic kinetic term, and leaves the Pauli-Feriz mass-term, £ {2} D C FF = h 2 — h^h^. The parameter choice 
X = — Y is called the Pauli-Feirz tuning, and will render massive gravitons ghost-free at linearized order on Minkowski 
backgrounds. In this case Ay^ = 0, -Byv = — X, Cyy = \X, Dyj = X, Eyj = — \X which is a special case of the more 
general situation discussed earlier. 

Rather than retain Lorentz invariance and be forced to use the Pauli-Feirz tuning to remove the ghost, it has been 
suggest that one can just fix the field = which implies that x = 0j removing it as a physical degree of freedom. 
Of course, this is not really a solution to the problem of the ghost, since we have just set the field to zero. However, 
as we will see in the next section that it is possible to impose a symmetry which is equivalent to this. The condition 
M^C/x = imposes an interesting structure upon the fields in theory when we use the transverse-longitudinal split 



language. Using this and (3.3 1 implies that k = —u^C,^ which can be differentiated to yield k = —u^Cu- This shows us 



that the K-field (i.e. the longitudinal component of the £ M -field) does not propagate. In stead , the n th time derivative 
of k is replaced by the (n — l) th time derivative of £ M . Evaluating the "kinetic scalars" (2.13) for the scalar k, yields 



(3.5a) 



(3.5b) 



Using (3.5b), the previously offensive term in (3.4) becomes 

(x + y)(Dk) 2 = (x + r)K^4c. + 2M Q CaV^v Al K + v^v AI KV Q v aK ), 



(3.6) 



and we observe that the multiple derivatives of k that are present are entirely spatial. The upshot is that there are 
no time-derivatives of the scalar k left, and, crucially no second time-derivatives of k in Ok. The term (X + Y)(\3k) 2 
in (3.4) is now longer problematic, and does not require removal. 
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IV. IMPOSING REPARAMETERIZATION INVARIANCE 



A key aspect of the theories under consideration here is the spontaneous violation of reparameterization invariance. 
It is interesting to see under what conditions it can be be reimposed on the theory. Therefore, we consider how 
the vector held £ p sources the perturbed gravitational held equations, and under what circumstances its components 
decouples from the field equations. From (2.1) the field equations for the metric are S^G^ = % , kGS e T iiv + SJJ^ , 



where ^ E J7 A " / is the dark energy momentum tensor and contains contributions to the field equations from the dark 
sector. In [5] we showed that 



(4.1) 



It is useful to note the terms in S E U fil ' which are present due to the background dark energy- momentum tensor U^ v . 
The components of S E U 11V are written as perturbed fluid variables, 



£ E t/% = Spu»u v + 2(p + P)v^u v) + SP^ V + PW V , 



where one can obtain 



p + U (2p - A w + 3B W ) 



X -(A W + p)x +(p + B w ) [ -h + diU 1 



(4.2) 



(4.3a) 



SP 



P + H(2P - B w + 3D W + 2E W ) 



(B w -P) X -l(P + 3D W + 2E W ) -h + diu 



(4.3b) 



(p + py = (p + c w )d i x + (p- c w )w i 



(4.3c) 



pn l j = 2{p - e- 



w) 



^h i i + d<tu j) -\& j (\h + d k J c ) 



(4.3d) 



These effective fluid variables define how the components of t he ve ctor field £ M sources the gravitational field equations. 
If one, or both, of the fields x an d u % does not appear in (4.3) then that field is not dynamical and hence can be 



completely ignored. It is clear that particular choices of the free functions in the mass-matrix it will be possible to 
achieve this. When one or both does not appear, it means that the theory is invariant under the symmetry associated 
with that field. Therefore, we can impose reparameterization invariance in three natural ways: 

• ^'-field decouples from the system when p + 3H(p + P) = 0, P + H{P + 3D W + 2E W ) = 0, A w = -p, B w = P, 
Cyy = ~P-, P = -Bwi P = CyVi P = -EWj D\v = ~ P- Hence, in the "fully" reparameterization invariant case, 
where the theory is forced to be invariant under — > £ M the only values of p, P that are allowed are those 
which are provided by a cosmological constant, p = —P, and all perturbed fluid variables vanish. Neither the 
X~ nor the a/-field propagate. 

• X — u ^ fi e ld decouples from the system when the parameters satisfy Ayj = —p, Bw = P, Cw = —P, 
p + H (2p - A w + 3B W ) = 0, P + H (2P - B w + 3D W + 2E W ) = from which we can deduce that p + 3H{p + 
P) = and P + T~L(P + 3_Dyv + 2P>w) — 0. These equations appear to leave two coefficients, Dyy and Ew, 
unspecified. If we now define two parameters /3 and p via -Dyy = P — P — \p and Ew = p + P then we find that 
(3 = (p+P)^f- which is the definition of the relativistic bulk modulus and p can then be interpreted as a rigidity 
modulus of an elastic medium. Hence, in the case where we impose time translation invariance, t — > t + x- but 



not spatial translation invariance, then we find that the theory must be EDE. The equations of motion (2.19) 
become 



-3H[P + 3f3H] X = 0, 



(4.4a) 



[p + Pp' l + [P + H{P + p)]uj l - [Ew - P]V 2 uj 1 - [E w + D w ]d l d k u k 
= -[P- Ew] dj0 + \[P + Dw] d l h. 



(4.4b) 
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Note that (4.4a I vanishe s for arbitrary values of x since = (p + P)P/p, and that there is a propagating vector 



degree of freedom, d 1 . ( 4.4b I is the equation of motion presented in 60 . In this case the mass-term for the 
gravitons is 



'{2} 



(w 2 ~lp)h 2 +2{w + p)h llv W v 



(4.5) 



where w = P/p,p = p/p. This case is equivalent to setting Aw = in the Minkowski space case. Since x does 
not appear in (4.3) it is no longer a physical degree and there is no ghost. 



• w ! decouples when p + B w = 0, P + 3D W + 2E W = 0, p = Cyy, P = E w from which we can deduce that 
Byv = — Cyy = — p and Dyj — —Ew = —P. Therefore, we see that in the case where we impose spatial 
translation invariance, x % — > x 1 + , but not time translation invariance, then the perturbations have some of 



the characteristics of massive scalar field theory as explained in [2]. The equations of motion (2.19) become 



[A w + p] X + [Aw + U(AAw + P - 3P)] * + [p + P] V 2 X 

+ [(A w + 3( A w - p - 2P)H)H + (A w + 4p + 3P)H] x = -\[p + P]h, 



(4.6a) 



[p + 3H(p + P)]d l X = 0. 



(4.6b) 



the entropy 



Note that (|4.6b|) vanishes for arbitrary values of x due to the background conservation equation. If we define 

'SP 



wT = 



Sp 



-w)5, 



and set Aw = [w + e(l + w)]p we find that 

r , 1 

WL = 



1 + e 



w(l + e) — 1 ' 



(4.7) 



(4.8) 



This implies that this theory is a scalar field theory with a non-standard kinetic term. The mass-term is given 
by 



£{2} — 



iwp 



h 2 - 2h^h^ 



(4.9) 



which we note does not satisfy the Pauli- Feirz tuning. This case is equivalent to setting Cyy = in the Minkowski 
space case. Since u/ does not appear in (4.3) it is no longer a physical degree and there is no ghost. 



V. DISCUSSION 



It is well established in the literature, that can be ghosts in general massive gravity theories, and various methods 
have been devised to remove them. If one imposes Lorentz invariance, then one is forced to use the Pauli-Feirz tuning 
to excise the ghost. If one is willing to give up Lorentz invariance, then certain parameter choices allow for ghost-free 
massive gravity theories. We have shown that the theory (2.2) with the (3+l)-decomposition of the mass-matrix 



(2.15) imposed with just time-translation invariance constitutes a linearized theory with healthy Lorentz- violating 



massive gravitons with 5 physical degrees of freedom. That theory is exactly the EDE model previously discussed in 
the literature |60j . In addition, if one imposes spatial translation invariance then there is another ghost-free massive 
gravity theory with 3 degrees of freedom. 

In terms of the EDE parameters and graviton masses, Aw = rn-Q — —p, Bw — —2m 2 — P, Cw — 4m 2 = — P, Dw = 
Ew — —2m\ = p + P. The masses can be conveniently parameterized by some overall mass scale 

Plp.fi P has units of mass squared), 



4m§ =P-P- 



M z = p> 



Hq (since we wrote = pu^u^ 



-M 2 



■\wM 2 , 



ml = —h(fi + w)M 2 



\(w 2 -Ip)M 2 



-\wM 2 . 



(5.1) 



We defined p, = p/p in analogy with w = P/p. 




Figure 1: The shaded region denotes the range of values of the equation of state w and shear modulus jX which yield sound 
speeds less than unity, which is where the inequalities ( |5.3[ ) are satisfied. The red (solid) lines denote lines of constant 
m\ = \ M 2 (0.6, 0.3, 0) from left to right and the blue (dashed) lines of constant m| = |A4" 2 (0.6, 0.2, — 0.2), again from left to 
right. 



To connect to dark energy, we note that the fraction of the total energy density that is dark energy is linked to 
the mass scale M 2 via fl de = M 2 /(3Hq). Hence, we see that the "natural" scale for the masses in order for the 
modification of gravity to act as a source of cosmic acceleration is of the order the Hubble parameter, and are all 
multiplied by order unity "corrections" denned by two parameters which encode the properties of the elastic medium: 
its equation of state parameter w and shear modulus fi. The longitudinal and transverse sound speeds of EDE are 

m 

<*=*> + d = T^- (5-2) 
1 + w 1 + w 

Stability and subluminality require that < c 2 < 1, so that we have the following constraints on the possible values 
of fj, 

- \w{l + w) < ft, < §(l-w 2 ), 0<fi<l + w. (5.3) 



In Figure [T]we plot the allowed values of (w,ft) which satisfy (5.3), and some lines of constant m 2 an d m\. Observa- 



tionally the values of w, M 2 and ft (only ft is the "new" parameter) can be constrained and then (5.1 ) used to obtain 
the graviton masses. 

Elastic dark energy and massive gravity share two common features. First, they are both constructed from rank-4 
tensors, (the elasticity tensor and mass-matrix, respectively) and these tensors have identical symmetries in their 
indices. Secondly, they both have five propagating degrees of freedom. The extra degrees of freedom in elastic 
dark energy may have a different fundamental origin to those in massive gravity, but they enable an interesting 
interpretation to be extracted from massive gravities. Our interpretation is that massive gravity is the manifestation 
of rigidity of spacetime. 
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